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or Dirichlet’s limit. Butif «and & are not identical, and 
if we know their difference with sufficient precision—that is, 
within a negligible fraction of the wave-length 2c” where 
the lower bound of m is as stated above—nothing is easier 
than to calculate A, and Ag, and thus F,(v). It will be seen 
that F(z), as æ increases through £, passes quite continu- 
ously from f(E—0) to f(E+ 0), since A, changes from 1 to 0, 
and A, from 0 to 1, in a perfectly continuous manner. The 
Fourier integral F(x) is seen to be continuous however 
discontinuous f(#) may be. F and f cannot therefore be 
identical. They are numerically equal in the sense that 
F —f is infinitesimal compared with f, and that too for every 
value of æ within the specified range of æ. It is possible 
that many ordinary mathematical theorems state only 
numerical equalities in the above sense, because in the proof 
of such theorems it is usual to neglect an infinitesimal 
at some stage of the argument. 


XXVII. On Simultaneous Operational Calculus. 
By BALTE. VAN DER Pou, D.Sc., and K. F. Niessen, D.Se.* 


L CCORDING to J. R. Carson + Heaviside’s operational 

L calculus can be based upon the transformation of 
a given function h(2) into its “image” f(p), related to the 
former one by means of 


f(p)=p\ ehle) de. 
Jo 
This relation between image and original we write 


SPER). 


If the image f(p) is known, one can obtain the original 
function A(s) in the way pointed out by Bromwich f, 
leading to one single function h(#) only on account of 
a theorem of March t. 

The number p, whose real part must be positive 


Rip) >0, 


* Communicated by the Authors. 

+ J. R. Carson, ‘Electric Circuit Theory and the Operational 
Calculus’ (McGraw-Hill, New York, 1926). 

t For references about these articles see Balth, van der Pol, Phil. 
Mag. viii. p. 861 (1929), which will be quoted as I. See also 
H. M. Macdonald, Math. Soc. of London, xxxv. p. 428 (1902). 
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is called the parameter which operationally belongs to z, 
although really 7 (not p) and æ are related to one another 


as image and original : 
I 


p 
We now wish to introduce a simultaneous operational calculus, 
where two variables æ, y are treated operationally by means 
of two parameters p, q respectively. 


la ) (fi(p)=li(@), 

pr | or in general ; 

Ai] | , 

ae \ fo(q)=ho(y), 
based upon 


Si(p) =p ( e-h (x) da, 
vo 


ZORN e- Whe(y) dy. 
vat) 

As an example we will determine the original function 
H(z, y) corresponding to the image: 


1 : 


waa H(#, y). 


Regarding for an instance g as a constant, we first investigate 
the interpretation h(«) of 
1 


Sul?) = Gap 


Applying a known theorem (see Carson or van der Pol, 
Le. (8)), we have 


h(x). 


Ae- hla), 


where c is any constant. 
Choosing g for c, we get 


p č la 
— s hlt) . . . . . (1 
p-q p~ 1e) (1) 


where the sign = is related to the transformation LM 
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We might better write 


A(w)=hfa), . . . . . +» 


to indicate that there is still the parameter g left. Now we 
have, keeping q still constant, 


a 
——Se™ Fw. G 
Pare (3) 
In the known theorem, that when 


fe) = hz), 
1 a 
= f(p) =f hajde, . . . .. (4 


it is tacitly understood that 


it follows that 


Í h(a) da << e”, 
0 


This is seen by performing the partial integration 


AN em M (£) dx= -{ h h(a) dæ d. e7? 


sgt ii hay dwt riip: 


Therefore the application of (4) to (3) is only permitted in 
the case 

a P> 
giving 


See a (i (dayr 
0 


` This multiple integral can be transformed into a single one 
in a known way, so that 


1. (7 ype- à 
oat oo ae DT CDTi... 


* This transformation follows also directly when the left-hand side- 
of (5) is written as 
? 1 1 


p—9 pr-l" p 


and is considered as a “ produit de composition.” 
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Comparing (5) and (1, a we find 


Writing ee, 
we obtain aac aapil, e-%g"-1 ds, 


an incomplete ny hs 

In order to find the required function H(z, y) we have to 
transform hg(x) as a function of q operationally into a 
function of y: 


a= Gay [ewes hen... © 


With reference to J, form. (9) we have 


0, for y<s, \ 
e= Boy) = >0. . 7) 
EAy)= ie for y>s. j 1) 


With the function F,(y), defined on the right-hand side 
of (7), we have from (6) 


He Y= capt Bets . . & 


If we wish to know H(#, y) for a set of variables 2, y, where 
a<y, we have in F,(y) of (8) always s<a<y, and therefore 
from (7) 

F,(y)=1 for s<a<y, 


so that | 
H(z, y)= E ds = a for <y. 
Seeking the value of H (æ, y) in the case y< w, we have 
(y)=1 for s<y, 
F,(y) =0 for y<s<a, 
so that (6) becomes here 


Hl )= Gay) om ds=*, for y<a. 


We can state the result in the following way :— 


1 | for a<y, | 
Gree He N=) a ‘ (9) 
(al for ysa. j 
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For n=1, H(æ, y) would become a kind of a “roof-function” 
when we measure H(z, y) along the ordinate < above the 
region for positive values of x and y (see figure). 

This geometrical interpretation is only given as an illus- 
tration, but will not be used more in what follows, where 


we wish again to draw attention entirely to the analytical 
side of the problem. 


2. We will now consider an application of the simultaneous 
operational calculus to the evaluation of some integrals. 
Suppose we require the value of 


sl PE 2 < aw « 10) 


bps 


Writing 7 
u=23VĒ 2. .... AD 


we wish to know 
2s { I2 VEI VE gal. 
0 


The reason for this substitution (11) is to bring the Bessel 
function and its argument in a form suitable for operational 


interpretation. In I. form. (29), namely, it is proved that 
when p and « belong together in 


1 
ts 


p 
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the following operational transformation exists : 


n 1 
aiJ (2e) =p e P. a n  . (12) 
From the above it will be clear that we better write 


° EJ, (2V E) ETIVE) 
28= ( geri dé. 
Now for what follows the essential point is that we do not 
transform a function of & into a function of p, but that we 
first introduce new variables « and 8, not instead of E, but 
entirely apart from §, namely, by constructing the function 


Sa, s= Y (a£)? J (2V ak) (8E) 2Ja(2/ BE) dé. (13) 


ett 
For «=B8=1 we have 
31, t =28. . . . . . . (14) 
To evaluate S,,g as a function of'a and 8 we will use our 
operational methods : 


a function of æ will be transformed into a function of p, 
-L 


gr 3 
and a function of 8 into another of q, 


all 
B= 7 
From the mentioned theorem (see (12), where p and æ belong 
together) we know 
1 


atJ (2V a) pute Pt. n n (15) 


BEIKLVB grei. o. o. (16) 


In these operational relations & can be treated as a constant. 
Now, using for an instance the notation of I. form. (5), 


it follows from 
Fip) =ke) 
that 


F(Z) = lsa), 


if s is a positive number. 
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Applying this with s= and «=a or 8, we find’ from 
(15, 16) 
2 mi š p = _é 
CAPIVA = (F) Tr, . . (15a) 


s£ 
q 


(BB ?I(2 VBE) = (4) "e (16a) 


Therefore from (13) 


DATORE 


grt: dE 


1 gu ; 
p deed 
= aa), OT Dee a 


The symbol = contains here the aoe as well as the 


-== 8 tranformation. 


Substituting 
1 Wy 
ptg F=% 
we have 
2g L 9 M(n—1) 
ens s"—lds = 
eo, (w +g)" 


For the evaluation of Sag we need, therefore, the inter- 


: 1 : . 
pretation of pro” which was found in the foregoing 


paragraph. 
Applying (9), we obtain for Sa, g 
(Mes) y if «<p | 
Z | EON i 
Sag = 4 kea (dT) 
i naD g, if B< | 
U Ila) z a) 
and 
1 
S1= n° 


On account of (10) and (14) we find, therefore, for the 
required value of the integral S: 
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Meanwhile we found even more than that, since from (13) 
and (17) it follows that 


n 


TAR if a<8ß, 


v0 


i 
(° DOVME rE | 
| 


If we substitute i? 


aa, B=b, 
we obtain a special case of the discontinuous integral of 
Weber and Schafheitlin (comp. Watson, ‘ Bessel Functions,’ 


§ 13, 42), namely 
l/a « 
DII ip E E) iË a<b, 
(OD (aC) 
a0 t Bi 1 b n : 
| zala) > if b<a. 


3. Another application of simultaneous operational calculus 
can be made in the proof of the orthogonality relation for 
Bessel functions of the same order. 

As is known, this relation is 


‘1 
í uJn(au)Jn(bujdu=0, 
0 


where a and b are two different roots of 
Jaz) =0. 


This has been found from the following equation :— 
1 
(a?—b2) f uJ,(au)J„(bu)du = a3 y()Tn+1(4) —bT (a) Ty (0). 
. . (18) 


This relation can be proved by operational methods as 
follows. Substituting 


a=,/ a, b=vV P, ux 2V E, 


we have to prove 


-of CE 2J (2V aE) (BE)? Ja (21/ BE) dé 


n 
n+l 


a ia las n n nt 
= ym ” Jay (V B Ja B) — basJ, (a8 E dan 
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If we relate 


= >») 


g 
and use again the equations (15 a) and (16'a), the right- 
hand side of (19) is transformed into 


1 
ro R 1 1 
mai) ooo 0) 


For the transformation of the left-hand side of (19) we 
need the rule mentioned in I. form. (12), namely, it 


follows from 
F) = h(x), 
that 


(-5) 10) = aw" h(a), n>0. 


This gives in our case (n=1, c=a or £): 


ae cove 2-9 8 (ES?) 


senin eV 5B) =~ a2 (PS). 


so that the left-hand side of (19) becomes 


ee 
Pel = ee 
(21) 


The equality between (20) and (21) is directly proved by 
introducing 
1 1 
s=- +) 
P 9 
and hence the required proof of (18) is given in an 
operational way. 
Natuurkundig Laboratorium der 


N. V. Philips’ Gloeilampenfabrieken, 
Eindhoven, July 1930. 


